arXiv:1503.08887v2 [cond-mat.str-el] 21Jul2015 


Theory of excitations and dielectric response at a spin-orbital quantum critical point 


Daniel Ish 

Physics Department, University of California, Santa Barbara, California 93106, USA 

Leon Balents 

Kavli Institute for Theoretical Physics, University of California, Santa Barbara, California 93106, USA 

Despite possessing a local spin 2 moment on the iron site and a Curie-Weiss temperature of 45 K, 
the A site spinel FeSc 2 S 4 does not magnetically order down to SOrnK. 1 2 Previous theoretical work 3 ® 
by Chen and Balents advanced an explanation for this observation in the form of the “Ji- A” model 
which places FeSc 2 S 4 close to a quantum critical point on the disordered side of a quantum phase 
transition between a Neel ordered phase and a “Spin-Orbital Liquid” in which spins and orbitals are 
entangled, quenching the magnetization. We present new theoretical studies of the optical properties 
of the J 2 -A model, including a computation of the dispersion relation for the quasiparticle excitations 
and the form of the collective response to electric field. We argue that the latter directly probes a 
low energy excitation continuum characteristic of quantum criticality, and that our results reinforce 
the consistency of this model with experiment. 


I. INTRODUCTION 

A suite of experimental probesPSEHi identifies the A-site 
spinel FeSc 2 S 4 as a rare example of a orbitally degener¬ 
ate antiferromagnet which resists magnetic or orbital order 
down to a temperature of tens of millikelvin, ma kin g it a 
truly quantum paramagnet. It has been suggested^® to lie 
close to a quantum critical point, making it a compelling 
object of study. This prior theoretical work proposed the 
“J 2 -A” model for this compound in terms of spin two, S j, 
and spin one half, T ; , operators on the diamond lattice 
with the Hamiltonian 

h = j 2 y j s 1 -Sj + J2^~ b Y,s; a ) 

<M> * * 


where (i,j) stands for next nearest neighbor bonds and the 
on-site Hamiltonian TC/ is given by 


W = -- 


(V3 T* [(S?) 2 - (Sf) 2 ] + T* [3 (S’) 2 - Sf] ) 


( 2 ) 


with J 2 > 0 and A > 0. We have included in Eq. 0 an 
external magnetic field B , taken for concreteness along the 
crystalline (001) axis, whose effects we will study further in 
the following. One should interpret this Hamiltonian as de¬ 
scribing the low-energy dynamics of the 6 d electrons on the 
Fe 2+ sites. Due to the tetrahedral crystal field, the d man¬ 
ifold splits into a lower e doublet and an upper t 2 triplet. 
These states are then filled in a high-spin configuration on 
the assumption that on site Hund’s Rule exchange domi¬ 
nates the crystal field splitting, giving an overall spin 2 (Sj) 
together with a two fold orbital degeneracy (T*). The J 2 
term in the Hamiltonian is a NNN antiferromagnetic ex¬ 
change term which arises in the standard way from virtual 
hopping between Fe 2+ sites and the A term represents the 
effect of spin-orbit coupli ng at second order, coupling the e 
hole to the overall spin 2 p®M 

With J 2 = B = 0, 5C° describes a system of uncoupled 
Fe 2+ sites with tetrahedral geometry under the influence 


of spin orbit coupling. This splits the 10-fold degenerate 
high-spin manifold into five equally spaced levels separated 
by energy A. These are, in order of ascending energy, an 
ai singlet, a triplet, an e doublet, a t 2 triplet and an a 2 
singlet. The ai ground state takes the form 

_L j x 2 _ y2 ) | 0) + 1 |3,2 _ r 2) ( | 2) + |_ 2)) (3) 

where the number in the second ket in each product refers 
to the S z eigenvalue. Critically, this state has zero aver¬ 
age magnetization along all axes. In the presence of J 2 , 
single site ti excitations (“triplons” or “spin-orbitons”) ac¬ 
quire a k-dependent dispersion, but remain massive until 
J 2 /X = 1/16 at which point the system undergoes a quan¬ 
tum phase transition to antiferromagnetic ordering at w ave 
vector q = (27 t, 0, 0) and symmetry related wave vectors.^® 
In the J 2 -X model, this ordering actually happens indepen¬ 
dently on each of the fee sublattices. A NN exchange term, 
Ji, is also all owed by symmetry,®^ but is expected to be 
much smallei®^ 12 and is difficult to distinguish from J 2 ex¬ 
perimentally at k = 0, which will be the regime of focus in 
this article. This term controls the relative orientation of 
the magnetizations of the A and B sublattices. 

One can argue for the consistency of this model with 
FeSc 2 S 4 at lowest order in terms of two distinctive experi¬ 
mental observations: the lack of observed magnetic order¬ 
ing down to temperatures of 50 mK® 2 ] and the observation 
of a low energy mode at momentum (27r,0,0) by neutron 
scattering.® Taken together, these suggest that the model at 
least qualitatively predicts the behavior of the compound 
if we suppose that the compound lies on the disordered 
side of the transition close to the QCP. Indeed, previous 
estimates 3 ® of the magnitude of these couplings in FeSc 2 S 4 
give A ~ 22.1 K and J 2 ~ 1.29 K, giving J 2 /X ~ 1/17 and 
putting the material just on the disordered side of the tran¬ 
sition. 

Given the striking properties of this material, it is no sur¬ 
prise that it has been the subject of recent experimental in¬ 
vestigations, primarily focused on its optical properties.®^ 
In this article, we attempt to make contact between this 
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FIG. 1: Spectrum of with magentic dipole allowed 
transitions (blue) and electric dipole allowed transitions 
(red, dashed). The number of lines in a level indicates the 
degeneracy of that level. 


model and the observed optical properties of FeSc 2 S 4 in 
order to argue for its continued consistency with experi¬ 
mental results. In order to do this, we will both investigate 
the fate of single site dipole allowed transitions in the pres¬ 
ence of non-zero J 2 , predicting both the shift in location 
and g factor, and investigate the character of the collective 
response to k = 0 electric fields. We also derive expres¬ 
sions for as-yet unobserved quantities and propose possible 
experiments to measure them. 


II. MAGNETIC DIPOLE EXITATIONS 

The single site problem has a single magnetic dipole al¬ 
lowed transition from the ground state, a\ — > t\. We expect 
that in the presence of J 2 , this excitation will still result 
in a peak in the AC response of the material to magnetic 
fields, but will shift in energy. Given the proposed prox¬ 
imity of the material to a quantum phase transition, we 
expect perturbation theory to be inaccurate when predict¬ 
ing the location of these excitations. To see this, note that 
corrections to the energies of low lying excitations due to 
J 2 must be comparable to A close to the QCP, since the 
J 2 = 0 gap (of magnitude A) has almost closed. Thus, to 
determine their dispersion, we instead calculate the RPA 
susceptibility in the presence of a (001) directed field and 
investigate its pole structure as a function of k. 

As we will see in the present section, this corresponds to 
treating the t\ excitations as non-interacting boson excita¬ 
tions (triplons) above a variational ground state (the mean 
field). On a practical level, this means that we will keep a 
smattering of selected terms at higher orders in J 2 . This 
approximation ultimately relies on the fluctuations about 
the mean field being small, i.e. on the triplons being di¬ 
lute. This approach, of course, has its own dangers. For 
one, the astute reader will notice that this approximation 
also breaks down as we approach the critical point, where 
these triplons condense. The hope here is not that the RPA 


is perfect close to the critical point, simply that it does bet¬ 
ter than perturbation theory due to the inclusion of these 
higher order terms. The natural second concern is that our 
choice of a description in terms of free bosons was incorrect, 
i.e. that we have chosen an unfortunate auxiliary field. Put 
another way, we may have decided on a bad prescription for 
which terms of the perturbation series to include. As some 
modest check on this, we will compare the RPA results with 
those of low order perturbation theory and consider more 
carefully the validity of this description wherever the two 
disagree. 

A. Formalism and B = 0 Magnetic Dipole Excitations 

We consider the imaginary time dynamic susceptibility 
Xij (n - r 2 ) = (T T Sf (ri)S^(r 2 )> (4) 

where 

§"(r) = 5 "(t) - (S 1 *) (5) 

Performing a Hubbard-Stratonovich transformation to de¬ 
couple the exchange term, we find the partition function to 
be given by 

Z = j D[$]e _Se// ^ (6) 

with the effective action for the auxiliary field 

Sef f m = \j dT^J^i ■ Tj - In W[*] (7) 

where IF[<f>] is the partition function for the J 2 = 0 problem 
with a magnetic field of —applied to each site (the 
“single site” problem). The dynamic susceptibility is then 
related to the propagator for the auxiliary field, fi(k, w), 

by 

x(k,w) = J(k) -1 — J(k) _2 fi(k, w) (8) 

where J(k) is the Fourier transform of the interaction 

J(k) = J 2 cos ( k ■ A ) (9) 

A 

with A the 12 fee nearest neighbors. If we then expand the 
action to second order in $ about its saddle point (d> = 0), 
we find that the bare propagator for the auxiliary field, Oo, 
is given by 

Ho = [J ( k ) -1 + X'oM’" ' (10) 

with Xo(w) the dynamic susceptibility of the single site 
problem and the inverse being, of course, the matrix in¬ 
verse. Together with Equation [8] this gives the RPA sus¬ 
ceptibility 

Xrpa{ k,w) = XoMt 1 + ^(k)xoMr 1 (11) 
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which we analytically continue to extract the real-time RPA 
susceptibility in terms of the real-time single site suscepti¬ 
bility. We can write the real time single site susceptibility 
in terms of the spectral representation 


xrM =Y 1 

o^o 


(0 \s»\j) <OTQ) 

Ej — E 0 — ui 


(om?) u\sm 

Ej — Eq + u) 


with |0) denoting the ground state. Using this, we find that 
both the single site and RPA susceptibilities are multiples 
of the identity and that the RPA susceptibility exhibits 
poles at 


w(k) = A^l + ^J(k) (13) 

which predicts a pole in the k = 0 susceptibility at 

w(0) = A^/l +48y « 1.95A (14) 

at the predicted value for J 2 , J 2 = A/17, as previously 
reportedP Since we have neglected all interactions, robbing 
the triplons of any decay channels, these poles lie on the 
real-frequency axis and give rise to S function peaks in the 
imaginary part of the real-time susceptibility. Expanding 
this dispersion to first order in J 2 gives 

w(k) = A + 2 J(k) + 0( J 2 2 ) (15) 

in agreement with first order perturbation theory. 


B. Magnetic Dipole Excitations with B + 0 

The preceding analysis can be repeated in the presence 
of a magnetic field with little change. The location of the $ 
saddle point simply shifts due to the presence of a magnetic 
term in the associated single site problem, leading to a 
different mean field. The new saddle point is of the form 

$ 0 (r) = -%p (16) 

with /3^ determined by the mean field consistency equation 

P = 12 J 2 (S») 0 (17) 

where (-) 0 stands for averages taken in the single site prob¬ 
lem in the presence of the field (B^ — with 

the applied field. In our case, since the applied field is con¬ 
sidered only along the (001) direction, we have that only 
/3 Z is non-zero. The form of the RPA susceptibility is the 
same as in Equation 0 save for the fact that the single 
site susceptibility is now calculated in the mean field 
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FIG. 2: Energy of the t\ triplet excitations versus B in 
the single site problem (dashed) and as given by RPA at 
k = 0 (solid). Field and energy are both measured in 
units of A. 


not be well captured by the the J 2 -X model and have not yet 
been experimentally observed. In the analysis that follows, 
the poles due to these transition exhibit some pathological 
behavior. In particular, as B —» 0, they return to the single 
site energy of the t 2 excitations, in contrast with the be¬ 
havior of the ti excitations. This is also in contrast to the 
perturbative result, which gives a shift of the t 2 energies at 
second order in J 2 . Consequently, we will not present the 
results of the B / 0 RPA for the two t 2 excitations that 
become allowed. 

For B ^ 0, the analysis of the pole structure of the RPA 
suceptibility is complicated significantly by that fact that 
the single site susceptibility is no longer a multiple of the 
identity. We only retain that the susceptibility is block- 
diagonal in the zz and xy blocks. Even with this reduced 
symmetry, there is still some remaining structure that we 
can exploit. Using the spectral representation, we see that 
Xq x = Xq V an d Xq V = — XcT ■ Together with the fact that 
X( X is purely imaginary, this implies that the single site 
susceptibility is diagonalized at all frequencies (and hence 
all imaginary times) by the same unitary transformation. 
The eigenvalues of the xy single site susceptibility are then 


X^(r) = X x 0 x (r) ± ix xy ( r ) = \ {T t S*(t)S ± ( 0)) q (19) 


and, of course, x zz • Together with Equation EH gives that 
the eigenvalues of the RPA susceptibility are 


Xrpa^ w ) 


Xq M 

1 +J(k)xo (w) 


( 20 ) 


= (T r ^(r)§^(0)) 0 


(18) 


Applying the magnetic field reduces the tetrahedral sym¬ 
metry of the single site problem, causing some transitions 
«i —> t 2 to become allowed. As will be discussed in Sec¬ 
tion III the shifts in the energies of the t 2 excitations may 


and 


Xrpa^i w ) 


xg*M 

i + j(k) x rM 


( 21 ) 


Furthermore, we can actually see by the spectral repre¬ 
sentation that Xo ( w ) = X'd( —1 w), which implies the same 
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FIG. 3: Energy of the t-\ triplet excitations versus k x for a 
few values of B, at J 2 = ^f. Field and energy are both 
measured in units of A. B = \ with A = 22.1 K 
corresponds to a physical fi eld of B v hvs = 8.23 T. (See 
Section III C 1|) . 


result for Xrpa • So, we actually need only investigate the 
pole structure of Xrpa an d Xrpa ■ This is quite easy to do 
for x 22 , since only one state has a non-vanishing S z matrix 
element with the ground state at all B. Calling this state 
| z) and its energy e z , we find poles of Xrpa 

w»(k) = e-y 1 + l(°l^ 2 l~>| 2 (22) 

where it should be noted that e z , |0), \z) and S 2 all depend 
on B through the mean field. This result can also be seen 
to be consistent with perturbation theory to first order in 
J 2 , though one must carefully track the dependence of the 
mean field on J 2 in order to obtain all the terms. 

The poles of Xrpa are more difficult to extract exactly, 
due to the larger number of states contributing to Xo ■ We 
can cast their location as the roots of an 8tlr order polyno¬ 
mial, however, and solve this polynomial numerically as a 
function of B and J 2 , producing Figures [2] and [3] Notice 
that non-zero J 2 actually reinforces the linear behavior of 
the triplet at k = 0 in small field. This is easily understood 
by noticing that the mean field is strictly smaller than the 
applied field, since the antiferromangetic interaction im¬ 
poses an energy cost to uniform magnetization. 


C. The Linear B Regime 

1. Computation of the g Factor 


where we have anticipated that this splitting may depend 
on wave vector, and inclu ded a factor of 2 to account for the 
fact that B = 2/j p B p h ys EUTD Now, if we define the total 
field felt by the site as 


B s = B - /3 2 


(24) 


we then find 


fl(k) 


o de+(B s , k) 
2 SB, 



dB s 

dB 


B =0 


(25) 


Implicitly differentiating the mean field consistency equa¬ 
tion gives 


dB s 

dB 


B =0 


1 

1 + 4 ^ 


(26) 


For the other derivative, we use first order regular per¬ 
turbation theory on the polynomial derived from RPA to 
find 


de+(B Sl k) 


dBs 


= 1+4 


J(k) 


(27) 


B s = 0 


so that 


S(k) 


A + 4 J(k) 
A + 4J(0) 


(28) 


At the zone center, this predicts no modification to the 
single site g factor of 1. Curiously, we also find that the 
g factor at the ordering wave vector decreases to zero as 
we approach the critical point. This agrees well with the 
qualitative behavior observed in Figure [3j Notice also 
that by B = 1/2 and J 2 /A = 1/17, we are already well 
outside of the the linear regime at the ordering wave vector. 


2. Small Fields 

In the linear B regime, we actually find that the system 
responds identically to a static applied field B = Bn in 
any direction. To see this, we show that the change in the 
ground state magnetization and dynamic susceptibility of 
the single site is isotropic to first order in B. Since these 
are the only two quantities from the single site problem 
that enter the calculation of the RPA susceptibility, this is 
sufficient to show that the response is isotropic at the RPA 
level. For this section only, B will refer to the magnitude of 
an arbitrarily directed field, rather than the magnitude of 
field applied along the (001) direction. So, let | in) denote 
any exact eigenstate of the single site Hamiltonian in the 
presence of B and write 


In order to characterize the splitting of the magnetic 
triplet in low field, we can compute an effective g factor 
in RPA 


ff(k) = 2 


de+(B,k) 

dB 


B =0 


(23) 


| m) = |too) + B |toi) + 0(B 2 ) 
For the ground state, we can see 

|0i> = jYl l^'°) Oo|n-S|0 0 ) 

0 


(29) 


(30) 
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giving 

m = = 2m [<°o|S|Oi>] = ^ + 0(£ 2 ) (31) 

(0|0) A 

so that the magnetization of the ground state is indeed 
isotropic to first order in B. 

As for the single site susceptibility, we will receive two 
corrections to first order in B , one from the first order cor¬ 
rection to the energies and one from the first order correc¬ 
tion to the states. We will write this as 


xS>> B) = x£>, 0 ) + B «» + <») + 0(B 2 ) 

(32) 

with r]E begin the first order correction from a shift in the 
energies and rj s being that from the states. Let us focus 
first on r] s . Investigating the relevant product of matrix 
elements from the spectral representation for y'o gives 


<0|S» (n|S ,! '|0) = (Oo|S>o) (nol^lOo) 

+ B(0 1 |5'*|n o ) (no|^|0o) 

+ B(0o|5Hni)(n o |^|0 o ) 

+ B(Oo|S'*|n 0 } (nil^lOo) 

+ B (OolS^Ino) (n 0 \S , '\0 1 )+O(B 2 ) 

(33) 


to first order in B. Since some expectation of the form 
(no|5“|0o) appears in each of the summands, we can see 
that this vanishes to first order in B for all states except 
those that evolve from members of t\. Furthermore, for 
\no) G ti, we see that |m) is orthogonal to t\. Since S a |0o) 
lies entirely in 1 1 , only the first, second and last summands 
contribute. Using Equation |30[ we find 


ry 


IJ.U _ 
S 


4 i 

Xa 2 



n a e 


OiflV 


(34) 


where n = n a e a and e is the Levi-Civita symbol. Indeed, 
this contribution to the dynamic susceptibility is isotropic. 

For rjE, the matrix elements in the spectral representa¬ 
tion are all between zeroth order eigenstates, so we can 
again restrict our attention to the t\ states. Here we diag¬ 
onalize the perturbation (B • S) restricted to t\ to obtain 
the zeroth order eigenstates and the first order energies. 
Only two states receive corrections to their energies at first 
order. A straightforward computation then produces 


LLU 

Ve = 


4* Xlo 

(A 2 - oj 2 ) 2 


n a e 


Ctfl V 


(35) 


which we can easily see is also isotropic. 


III. RESPONSE TO ELECTRIC FIELDS 
A. Electric Dipole Excitations of a Single Site 


dipole allowed transition from the ground state ai —> t 2 . 
One would imagine that the story for these excitations 
ought to be similar to that for the magnetic dipole exci¬ 
tations: the peaks in the permeability due to these excita¬ 
tions will persist in the presence of exchange with a shift in 
position. In this case, it is actually considerably more diffi¬ 
cult to make these statements quantitative for a number of 
reasons. Perhaps most glaringly in contrast to the question 
of magnetic dipole excitations, we do not know the identity 
of the operator which couples the single site problem to an 
electric field (P). One can demand that such an operator 
transform as a vector under the point group, i.e. as t 2 , but 
this still leaves the magnitude of its five reduced matrix 
elements undetermined. 

Furthermore, we have no reason to believe that the shift 
due to the J 2 term ought to be dominant over those due 
to all of the other symmetry allowed exchange terms ab¬ 
sent in the A model. Though J 2 is thought to domi¬ 
nate the other couplings, it only corrects the energy of the 
electr ic t riplet at second order. Other symmetry allowed 
termsP® e.g. (T?Sj) ■ (T^S,;), give corrections to the en¬ 
ergy at first order in their coupling constants that, taken 
together, might dominate those of J 2 . Since there has not 
yet been an experimentally unambiguous observation of the 
electric triplet excitation, we abandon the question of quan¬ 
titatively predicting the shift in its energy. 


B. Collective Response and Critical Behavior 

1. Coupling the Critical Theory to Electric Fields 

We can also consider the form of the response to electric 
field coming entirely from the low-lying magnetic excita¬ 
tions. As we will see, multiple triplon excitations possess 
the correct symmetry to be produced through coupling to 
the electric field. To determine the contribution of such 
processes to the electric field response of our model, we 
first restrict our considerations to low energy modes near 
the ordering wave vectors. Either expanding Ho close to the 
ordering wave vectors and for small frequencies or perform¬ 
ing a symmetry analysis® produces a Gaussian theory of 
the form 

Beff['lp a J=i / 7^ G M 1 ( k ’* W ™)|' , /’a, M ( k UW n )| 2 

^ a : fi : cu n ' ' 

(36) 


with 

G“ x (k,io;„) = -(iu n ) 2 + kV^k + r 2 = -(*w„) 2 + e 2 (37) 

where /i labels which ordering wave vector each of the fields 
came from, and a labels the sublattice. The matrix V), is 
of the form 

fv-t 0 0 \ 

14 = [ 0 V2 0 (38) 

\0 0 v 2 J 


Since the tetrahedral symmetry of the single site problem 
does not include inversion, there is actually a single electric 
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with V y and V z obtained by permutation. if> is an order pa¬ 
rameter for the staggered magnetization, though depending 
on how we obtained this theory, \j) a may not be precisely 
the staggered magnetization on the a sublattice at the /i 
ordering wave vector. For one thing, we have rescaled the 
field in order to set the coefficient of the w 2 term to 1. Addi¬ 
tionally, expectations of $ are not precisely those of S (c.f. 
Equation [ 8 ]). We expect that such a theory should describe 
our system correctly on energy scales small compared to 
the magnetic bandwidth. 

Now, rather than investigate the microscopic origins of 
the coupling of an electric field to this model, we simply 
investigate which couplings are allowed by symmetry. We 
expect couplings through a term linear in the applied field, 
E a P a , where P a is some function of the order parameter. 
This gives rise to a contribution to the k = 0 electric sus¬ 
ceptibility through the standard linear response formalism 

X? (0, w) = (P a (0, co)P b (0 , -w)) - ( P a ) < P b ) (39) 

We do not have couplings before second order in the or¬ 
der parameter, since the order parameter is odd under time 
reversal while electric fields are even. Beginning at second 
order without derivatives, we consider what restrictions re¬ 
quiring that a coupling of the form 

C abc E a ip b il) c = E a P£ (40) 

to transform trivially under the space group places on the 
tensor C, and hence the polarization Pq. We find that the 
symmetry allowed coupling is given by 

PS =ci ('.1.WA, 



(41) 


where ci, Ci and C 3 are undetermined by this analysis. The 
other components are related to the x component by simul¬ 
taneous permutations of the vector and wave vector indi- 
cies. Notice that Pq is odd under interchange of the sublat¬ 
tices, since inversion acts only to interchange the sublattices 
and gives a sign on the electric field. It can also be shown 
by expanding the effective action for $ given in Equation [7] 
that a coupling of the electric field to each site through the 
single site polarization operator (P) produces a coupling to 
the critical theory precisely of this form with c\ = c-i and 
c 3 = 0 . 

The next lowest order contribution to the electric field 
response should come from a coupling of the form 

D abcd E a %p h d c ip d = E a P1 (42) 

Due to the larger number of indicies, Pf contains many 
more terms than Pq so we omit a detailed discussion of its 
structure. 


2. Computation of the Electric Susceptibility 


Neglecting all couplings of higher order in fields and 
derivatives, we find the electric susceptibility is given by 


Xe = XeO + Xel (43) 

XeO = (PS( O, W )P 0 b (O,—cu)) (44) 

Xl\ = (P?( 0 , u,)P b ( 0 , -a,)) - (Pf) (P b ) (45) 

since (Pq) = 0 and 

<P 0 “(0,w)i^(0,-u;)> = 0 (46) 


since the internal momentum sum is odd under q —► —q 
while the Green’s function is even. Exploiting the particu¬ 
lar form of Po, we find that XeO is a multiple of the identity 
and 

Xeo(°’* w n) = ^ j 7^3 ^(q,*^ + Wn)Gx(-<l, ~^n) 

(47) 


with a = 2cf + 4c| + 4c§. Performing the Matsubara sum 
and analytically continuing to real frequencies gives 


xS(o 


dq 


coth 


(^) 


(27r) 3 e x (q) (4e x (q ) 2 - (w + iS) 2 ) 


Allowing S —> 0 produces 


Im(Xeo(°> w )) «sgn(w) 


d 3 q 

(2tt)3 


coth 






<KM - 2e x ) 

(49) 

(50) 


oc 0(|w| — 2 y/r) coth(/3w/4) 

Now, as for Xel , we can see on general grounds that 


ibcd 


Xel(° PUn) = 


P 


E 


_dq 

(2tt) 


(51) 


for some fantastically complicated tensor 77 . Terms with 
momentum dependence of the form q c qd for c 7 ^ d van¬ 
ish due to the momentum integration, together with the 
q —> — q symmetry of the Green’s function. The restriction 
that both Green’s functions come from the same ordering 
wave vector comes from demanding that Pi be invariant 
under the primitive lattice translations, together with the 


subtraction of (Pi). On identical grounds to Equation 49 
we then see 


Im(Xei(0,w)) oc sgn(w) 


(2tt) 3 a(q)^ 


<KM - 2 e b ) 
(52) 


a 0(|w| — 2 — 4 r) coth(/3w/4) (53) 
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FIG. 4: Dielectric loss due to Equation 50 as a function of 
wavenumber, in arbitrary units, at 5 K and 80 K. Both 
functions are normalized so that their value at 20 cm -1 is 
equal to 1. 


Arb. Units 



FIG. 5: Dielectric loss due to Equation |53| as a function of 
wavenumber, in arbitrary units, at 5 K and 80 K. Both 
functions are normalized so that their value at 20 cm -1 is 
equal to 1. 


IV. DISCUSSION 


incident light polarized both along the static applied field 
direction and transverse to it. At zero field, they observed a 
peak at ss 4.5meV in both polarization configurations. As 
field increased, the peak in the transverse direction split 
into two peaks in an approximately linear manner while 
the peak in the longitudinal direction remained unaffected. 
Fitting a line to the splitting of the peaks in the transverse 
susceptibility gave them g « 0.92. 

Following the discussion for Section |il B[ this is precisely 
what we would expect for the case of the field along the 
(001) direction. We saw that the longitudinal susceptibil¬ 
ity (i.e. \ zz ) received contributions from only one state 
and, expanding Equation |22| to first order in B , the energy 
of this state is independent of applied field to first order. 
The transverse susceptibility (i.e. the x-y block) received 
contributions from two states of the lower triplet, whose en¬ 
ergies split with a g factor just below the single site value 
of g = 1, due to the small wavevector of the incoming light. 
A priori, we might be surprised that the analysis with B 
along the (001) direction fits so well to data taken on a 
polycrystaline sample, but the results of Section |Tl C 2| tell 
us that this is exactly what we should expect, provided 
we are within the linear regime with respect to the static 
applied field. 

We would be remiss if we did not take this time to say a 
few words about the nature of the t\ excitations and what 
selection rules are relevant to this situation. It is tempt¬ 
ing to draw an analogy between the states of t± and those 
of a spin one triplet. Indeed, this analogy motivated the 
character of the analysis in Section |II B| After all, if we re¬ 
strict 0(3) to Td, the spin one representation becomes the 
t\ representation. So, at least formally, we can label the 
members of the triplet by m = 0 and m = ±1. In limited 
ways, this is even reflected in the response to applied field: 
the m = 0 and m = ±1 states are the zeroth order eige- 
nenstates with respect to the perturbation BS Z and the 
first order corrections to the energies are 0 and ±13/2, re¬ 
spectively. This is about where the analogy ends, however. 
The first order corrections to the eigenstates are non-zero 
(in sharp contrast with a true spin triplet) and there are 
higher order corrections to the energies of all members of 
the t\ triplet. 


A. Observations of the ai — > t\ Excitation 

1. Existing Observations 

Two recent THz spectroscopy experimentation FeSc 2 S 4 
have observed a well defined peak at in the range of 4.3meV 
to 4.5meV, in broad agreement with the previously esti¬ 
mated magnitude of A and Equation [14J which together 
predict that the t\ excitation at k = 0 should appear at 
approximately 3.7meVP21 Laurita et alP also performed 
this experiment in field and extracted results in remarkable 
agreement with the calculations performed here. To briefly 
recapitulate their story, they were capable of measuring 
the dynamic susceptibility in the presence of a field with 


Furthermore, one might be tempted to attempt to ex¬ 
trapolate selection rules from this analogy, saying that 
ai —> t\ is spin forbidden as a singlet to triplet transition. 
This would be incorrect, since in terms of the physical spins 
the transition is between different states of the S = 2 man¬ 
ifold, and no total spin change actually occurs. Since the 
electronic states see a reduced symmetry due to the crystal 
field and this reduced symmetry is communicated to the 
spins through spin orbit coupling, neither total spin nor to¬ 
tal orbital angular momentum nor total angular momentum 
are good quantum numbers on energy scales comparable to 
A, and we should not analyze selection rules in terms of 
them. The proper way to determine the selection rules for 
transitions between the eigenstates of TC/ is through using 
the Wigner-Eckart Theorem applied to Td. 
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2. Proposal for Future Measurements 


The authors are quite taken with the results of the com¬ 
putation of the g factor (Equation 28), and hope that it 
can be successfully measured soon away from k = 0. This 
would probably require an inelastic neutron scattering mea¬ 
surement on a single crystal. The benefits of such a mea¬ 
surement would be twofold. First of all, the result is un¬ 
usual and interesting in itself and it would be valuable to 
see it confirmed in the material. A successful fit of the g 
factor to the derived form would argue quite strongly for 
the predictive power of the J- 2 -X model for this compound. 
Secondly, and perhaps more importantly, a fit of the g fac¬ 
tor to Equation [28] would provide a measurement of the 
ratio J 2 /A, allowing one to estimate the proximity to the 
critical point directly. 


B. Continuum Weight in THz Absorption 


In addition to their observation of a transition matching 
the description of the a± —> t\ excitation, Mittelstadt et al. 5 
observed a curious continuum weight at low frequencies in 
the dielectric loss which was roughly linear at T = 80K 
and superlinear at T = 5K. One can attempt to explain 
this weight in term s of t riplon pair production using the 

focusing on the w 2 coth(/3w/4) 


IIIB 


analysis in Section 
term that appears in Equation |53[ At high temperature, 
coth(/?w/4) ss 4A 'bT/uj and we obtain the linear weight. At 
low temperature, the coth saturates and we obtain weight 
that grows as w 1 2 . Indeed, comparing Figure^ with the in¬ 
set in Figure 6 of Mittelstadt et al P, this functional form 
gives good qualitative agreement with the observed data. 
Within this picture, the fact that this continuum appears 
at low frequencies then becomes yet another signal of our 
proximity to the critical point, since we expect this weight 
only at frequencies in excess of twice the gap to triplon pro¬ 
duction at the ordering wave vector, as evidenced by the 
Heavyside 0 in Equation [53[ 


To be fair, there are a few objections that can be raised 
to this analysis. One could object that without a tem¬ 
perature dependent prefactor the w 2 * coth(/fu;/4) term fails 
to reproduce the observed temperature dependence of the 
absorption. The coth term decreases with increasing tem¬ 
perature, while the lowest frequency weight observed by 


Mittelstadt et al. increases with increasing temperature. 
The authors do not find this objection particularly com¬ 
pelling, as one would almost certainly find a temperature 
dependent prefactor upon a more careful analysis of the 
temperature dependence. Indeed, if one were to use a finite 
temperature version of the action in Equation [7] and expand 
about the ordering wave vector for small frequencies, one 
would find that the velocities in the Green’s function in 
Equation [37] depend on temperature. Repeating our analy¬ 
sis that led to Equation 53 but keeping more careful track 
of constants shows that the velocities appear among the 
( oj) constant prefactor that we have neglected. This is of 
course to say nothing of possible temperature dependence 
of the undetermined constants in the coupling described in 
Equation [42] 

Somewhat more concerning is the presence of the two 


( 1 const ) x coth(/3w/4) terms appearing in Equations 50 


and[53j which do not reproduce the observed uj dependence. 
The term in Equation [53] can be argued away consistently 
by claiming that the gap (y/r) is very small and we do not 
expect this term to be easily distinguishable in the presence 
of the large, dominating to 2 . The contribution from Equa¬ 
tion [50] cannot, to our knowledge, be so concretely claimed 
as negligible. One might hope that the presence of small 
but non-zero J\ might suppress this term by controlling the 
relative orientations on the A and B sublattices and lead¬ 
ing to some cancellation, but this turns out not to be the 
case. At the Gaussian level, J\ only acts to slightly reorga¬ 
nize our order parameters, shift the gap slightly and move 
the soft mode slightly away from the ordering wave vector. 
Pushing through the calculation, one finds only a change 
to the gap in Equation [50] It seems our only recourse is to 
argue that the coupling constants (ci, C2 and C3) ought to 
be small. While this is plausible, the microscopies of this 
coupling are quite daunting so we do not present a detailed 
analysis. 
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